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CHAPTER ONE 


FLOATING POINT INTERPRETER 


Introduction: 


The GRI Floating Point Intrepreter is a complete system that allows the 
user to process data in floating point arithmetic. Floating point arithmetic, 
through the we of multiple precision arithmetic and an exponential concept 
greatly extends the range of precision available to the user beyond that of 
fixed point arithmetic. It also, through utility routines, frees the user 
of the bookkeeping involved with scaling and unscaling of numbers that is 


necessary in a fixed point system. 


The GRI-905 has an instruction set which is known as machine language. 
The computer reads instruction words out of its memory and hardware is ac- 
tivated by the interpretation of each instruction word to cause the execution 
of that instruction. An interpretive software system fetches instructions 
which we shall call commands from the computer's memory and causes various 
subroutines to be entered as a result of the interpretation of the command. 
These commands fetched by the interpreter are also called psuedo-instructions 
because their format deviates from the machine's instruction format. The 


Standard machine format instruction is 


WORD 1 SDA MOD DDA 


WORD 2 [ADDRESS] (if a memory reference in- 
struction) 


A pseudo-instruction or command such as the ones used in the GRI Floating 
Point Interpreter looks like this: 


WORD 1 OP CODE 


WORD 2 [ADDRESS] (if a memory reference pseudo-instruction) 


The interpreter actually simulates the process used by the computer's 


hardware to execute an instruction. The interpreter fetches the OP CODE words 
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and addresses, sets up arguments, flags, and performs a function on the ar- 


gument(s) as specified by the OP CODE of the pseudo-instruction. 


An interpretive approach to floating point arithmetic provides the user 
with a functionally oriented language that makes usage of floating point 
arithmetic much easier than if it were done through a series of subroutines 
called in machine language. The user references floating point numbers 
with a single address which is the first address of the two word floating 
point number. The interpreter takes care of the address bookkeeping neces- 
sary for two word argument handling. The interpreter also maintains a set 
of accumulators much the same as an arithmetic unit. Arguments and results 
are manipulated and left in these accumulators. The interpreter utilizes 


two such accumulators plus an index register. 


There are a set of commands in the interpretive system that are not 
floating point arithmetic commands. These are program control commands such 
as conditional jumps and index register manipulators. The index is simply 
used to keep track of the number of times command loops are executed. These 
commands, although they could be effected by use of basic machine language, 
are also provided in the interpretive mode because they can save the user time 
that would be spent entering and leaving the interpretive mode, and almost 


always save space in terms of the coding needed. 


When the user is ready to execute commands in his program, he first is- 
Sues a machine language command that causes a jump to the interpreter to take 
place. The interpreter now assumes cancel and starts fetching commands which 
follow the jump that caused interpretive mode entry. If the user wishes to 
begin executing machine language instructions, he must issue an interpretive 


command that causes the interpreter to relinquish control. In essence, the 
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-machine is running in two different modes; a machine language mode and a 


psuedo-language mode -- in this case, a floating point language. 


The GRI interpretive system offers a novel error trap feature which 
may be invoked by the user to assist in tracking down places in the program 
where data values are causing error checks to occur. Errors such as dividing 
by 0, exceeding the capacity of the psuedo-accumulators in either the mantissa 
or exponent portions, etc., can all be caused by an unknown data base. A1l1l 
manipulations of data refer to manipulations in and out of the psuedo-accumulator 
called FAC. This accumulator behaves like the accumulator in an adding machine. 
It must be ieaded to initialize it, stored to save it, and all arithmetic op- 
erations leave their results in the accumulator. Commands with two operands 
are called binary commands and operate on a data word in user memory and 
the contents of FAC, replacing the result in FAC. Commands with one operand 
are called unary commands and operate on FAC, leaving their results in FAC. 


Let us consider a simple example: 


JU $SFI senter floating mode 

FLDA X sfetch X to FAC 

FMPY X 3X2 in FAC 

FSTA Tl 3store FAC in temporary loc 
FLDA Y 3fetch Y to FAC 

FMPY Y ;Y2 in FAC 

FADD Tl »x2 + Y* in FAC 

FSQT ; ‘[x2 + y2 in FAC 

FSTA R ;store result in R 


FEXT 3;exit from floating mode 
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1.2 Basic Package, SSFI: 


Floating point arithmetic capabilities are provided through an inter- 
pretive package. Associated with the package is an external<~> internal for- 
mat data conversion routine that can be easily tailored to the character set 


being processed. 


The interpretive package is invoked by a normal subroutine call. The 
call is followed by a string of commands that are established by use of 
equate statements during the assembly. The last command in the sequence 
causes a return to the calling program. Operations are performed using a 
pseudo acCumulator maintained locally by the interpretive package. The 
package also contains a 16 bit pseudo index to allow loops within the command 
sequence. Without this feature, it would be necessary to exit and re-enter the 
interpretive package and perform loop counts outside the interpreter, Al- 
though the latter procedure is, in most instances, faster in terms of time 
taken to do the loop, it usually involves considerably more code and, there- 


fore, takes more space. 


As an example of a typical problem programmed in the interpreter 
language, we evaluate the polynomial 
4 
Y = Ap + AX + ApX* + A3X? + AgX 
which can iteravely be expressed as Y = (( (A,X + A3)X + Ao)X + Ay)X + Ag 


as follows; 


LOOP: 


CONST: 


A3;: 


M4: 


JU 


FLDX 


FLDA 


FMPY 


FADDD 


FJIX 


FSTA 


WRD 


WRD 


SSFI 


CONST 
LOOP 


Y 


0,0 
A3-1 
X1,X2 
A41,A42 
A31,A32 
A21,A22 
Al11,A12 
A01,A02 


-4 
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;ENTER INTERPRETER 

;LOAD PSEUDO INDEX WITH ~4 

;LOAD PSEUDO-ACCUMULATOR 
;MULTIPLY IT BY X 

;DEFERRED ADD A3 (THEN A2, Al, AO) 
;COUNT THE LOOP 

;STORE RESULT IN Y 


3;EXIT THE INTERPRETER 


;STORAGE SPACE FOR ANSWER 
3;DEFERRED ADDRESS (GETS CHANGED) 
3TWO WORD FLOATING POINT VALUE OF X 


sFLOATING A4 VALUE 


3ONE WORD INDEX COUNT VALUE 
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1.3 Floating Point Format: 


Internal representation of a floating point number occupies two successive 
locations in memory and consists of a fixed point fraction (mantissa) with an 
associated exponent. The mantissa is in two's complement notation with a 
sign bit followed by 23 bits of significance. The binary point is assumed 
to be immediately to the right of the sign. The exponent, which is the power 


of two by which the mantissa is multiplied, has the range -200g to +1778 


(27128 to 2t127), this exponent is represented in "excess 200g 


notation by 
adding +200. to the true exponent. This requires a total of 8 bits and the 
range of the excess 200g notation is 000 to 377g, where 200g represents 20, 


Thus, a floating point number looks like: 


15 14 Bods eee 0 
Word 1 (16 bits) E | Fraction - most significant 
15 8 0. 
Word 2 (16 bits) | Fraction 
least significant Exponent _ 


This format allows an accuracy of 6+ decimal digits and a range of +1. 469 368x107 39 


to +1.701411x10+38, 


To obtain correct results, all floating point operations (except FLDA, 
FSTA and FNOR) require the floating point numbers being operated on to be 
normalized; that is, bit 14 of word 1 must be the most significant bit of 
the fraction (mantissa). The only exception to this requirement is a floating 
point zero, which has no significant bits--a normalized floating point zero 
is two words of all zero (mantissa = 0, excess 2008 exponent = 0). 
Note: The mantissa of a normalized floating point number other than zero 


has an absolute value in the range 1/24 mantissa; <1. 
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Examples: 

Decimal Internal Floating Point (octal) 
word 1 word 2 
1.0 | 040000 000201 
1.25 050000 000201 
-1.0 140000 000201 
-1.25 130000 000201 
100. 062000 000207 
-100. 116000 000207 
0.5 040000 000200 
0.25 040000 000177 
a 062207 166602 
TT /2 062207 166601 
-T 115570 011602 


1.4 Internal Registers: 


There are three pseudo-registers contained in the interpreter i) the 
pseudo-accumulator (FAC), ii) a temporary pseudo-accumulator (FTM), and 


iii) the pseudo-index register (FINDX). 


i) FAC - The floating pseudo-accumulator. This consists of three 
locations in the interpreter and is used to contain the left-hand argument 
of a binary floating point command as well as the results of any floating 
point command. It is organized as follows: 

FACHI - contains high order mantissa and sign of value in FAC 
FACLO - contains low order mantissa of value in FAC 


FACXP - contains excess 200g exponent of value in FAC 
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ii) FIM - temporary pseudo-accumulator. This consists of three lo- 
cations analogous to FAC. They are named FTMHI, FIMLO, and FTMXP. The 
temporary accumulator is used to hold an additional floating point value 
for those commands which require two floating point values in order to op- 


erate, e.g. a type II (binary) command (see 2.2.2). 


iii) FINDX - pseudo index. This consists of one location of the same 


name and holds the current value of the index. 


Note: FACHI, FACLO and FIMHI, FIMLO are treated as full 31 bit double 
precision quantities for the basic arithmetic operations add, sub- 


tract, multiply, and divide. 


2.1 Command Categories: 


FLOATING POINT MANUAL 
71-44-0001 


CHAPTER TWO 


BASIC COMMANDS 


The commands are of the following categories: 


Il 


II! 


eM 


VI 


load & store 


binary 


commands 


unary 


commands 


index 


commands 


conditionals 


exit 


> 


we 


we 


the command specifies the source or destina- 
tion of floating point data - the corresponding 
destination or source is the pseudo accumula- 


tor e 


the command specifies the source of the 
rightmost operand - the floating accumulator 
contains the leftmost operand. The result 


will be in the accumulator. 


the command merely specifies the function to 
be performed on the accumulator. The result 


Will be in the accumulator. 


the command specifies the source or destina- 
tion of an index value - the corresponding 


destination or source is the pseudo index. 


the command specifies an address to which con~- 
trol passes if the test defined by the command 
is true - the address must contain another 
floating point command. Tests may be per- 
formed on the floating accumulator, certain 


flags, and the index. 


this command causes a return to the calling 


program. 


262 
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The load & store (Type I) and binary (Type II) commands may 
specify deferred (indirect and auto-indexed) addressing mode. Deferred 
addressing in floating point commands operates exactly as in machine 


language. 


Command Descriptions: 


2.2.1 TYPE I COMMANDS -- LOAD & STORE COMMANDS 
LOAD FLOATING ACCUMULATOR (AC) 


mnemonic address code no. of words 


FLDA X Ol 2 


The contents of the location specified by X and X + 1 are treated as 
a floating point number and are loaded into the floating point pseudo ac- 
cumulator. The floating point number in locations X and X + 1 is split into 
three parts i) .X, which consists of the high order mantissa, goes into 
FACHI; ii) bits 8-15 of X + 1, which consists of the low order mantissa, 
goes into bits 8-15 of FACLO and bits O-7 of FACLO are set to zero; and 
iii) bits 0-7 of X + 1, which consists of the excess 200. exponent, goes 


into bits 0-7 of FACXP and bits 8-15 of FACXP is set to zero. 
DEFERRED LOAD FLOATING AC 


mnemonic address code no. of words 


FLDAD A 101 2 


The contents of location A is incremented by one, replaced in A, and 
the result is used as the effective address X; then the contents of A are 
incremented and replaced a second time forming the effective address X +l. 
The contents of X and X + 1 are then treated as a floating point number and 


loaded into FAC as explained under FLDA. 
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STORE FLOATING AC 


mnemonic address code no. of words 


FSTA X 02 Zz 


The contents of FAC are rounded into bit 8 of FACLO, bits 0-7 of FACLO 
are set to zero. Then FACHI, FACLO, and FACXP are packed into a floating 
point number and stored in X, and X +1. Note that this operation alters 


FAC so that it agrees with the value stored in X, and X +1. 


It is also possible for the rounding operation to cause exponent over- 
flow (excess 200g exponent exceeds +377.) . This can occur only if the 
number being rounded is very close to the largest possible positive float- 


ing point number. The value stored in this case will be X = 077777g, 


X+le= 177777.» and FXFLG will be set non-zero. A successful FSTA will set 


FXFLG to zero. 
DEFERRED STORE FLOATING AC 


mnemonic address code no. of words 


FSTAD A 102 2 


The contents of A are incremented twice as explained under FLDAD, form- 
ing effective addresses X and X + 1 into which FAC is stored as explained 


under FSTA. 
2.2.2 TYPE II COMMANDS —- BINARY COMMANDS 


All Type II commands depend on both FAC and the argument of the 
command to have normalized mantissas. If unnormalized numbers are: 
used, the results are unpredictable. A FNOR instruction (see 2.2.3) is 
provided to normalize any quantity if it is necessary to do so. Also, 


if all inputs are normalized, the results in FAC will be normalized as 
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will the value retrieved from FAC by use of an FSTA instruction. 


Type II commands can cause exponent underflow or overflow if the 
number created in FAC by the command has an excess 200, exponent out- 
side the range 0 to +377g respectively. The occurrence of either 
condition is indicated by FXFLG being non-zero after the operation 
has been completed. It may be tested by use of the FJEV command. 


The successful completion of a Type II command will set FXFLG to zero. 


FLOATING ADD 


mnemonic address code no. of words 


FADD . X 03 2 


The floating point number in locations X and X + 1 are added to 


the contents of FAC, and the result replaces FAC. 
DEFERRED FLOATING ADD 


mnemonic address code no. of words 


FADDD A 103 2 


The contents of A are incremented twice as explained under FLDAD, 
forming effective addresses X and X + 1, the contents of which are 


added to FAC, and the result replaces FAC. 
FLOATING SUBTRACT 
mnemonic address code no. of words 


FSUB X 04 2 


The floating point number in locations X and X + 1 are subtracted 
from the contents of FAC, and the result replaces FAC. 


DEFERRED FLOATING SUBTRACT 


mnemonic address code no. of words 


FSUBD A 104 2 
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Effective address is formed from A as in FADDD. 
FLOATING MULTIPLY 
mnemonic address code no. of words 


FMPY X 05 Z 


FAC is multiplied by the floating point number in X and X +1. 
The result replaces FAC. 


DEFERRED FLOATING MULTIPLY 


mnemonic address code no. of words 


FMPYD A 105 2 


Effective address is formed from A as in FADDD. 


mnemonic address code no. of words 


FDIV xX 06 2 


FAC is divided by the floating point number in X and X +1. The 
result replaces FAC. Divide check will occur if X, X + 1 is zero or 
not normalized. This causes FAC to be set to the largest possible 
floating point number of the sign which would be the result of the 
divide if it could take place, and the divide check flag (FDFLG) will 
be non-zero. A successful divide sets FDFLG to zero. 

Note - if both FAC and X are 0, the result will be the largest 
possible positive floating point number in FAC with FDFLG set non-zero. 

DEFERRED FLOATING DIVIDE 
mnemonic address code no. of words 
FDIVD A 106 2 


Effective address is formed from A as in FADDD. 
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FLOATING ADD MAGNITUDE 


mnemonic address code no. of words 


FADM X 07 Z 


The absolute magnitude of the floating point number in X and 


X + 1 is added to FAC. The result replaces FAC. 
DEFERRED FLOATING ADD MAGNITUDE 


mnemonic address code no. of words 


FADMD A 107 2 
Effective address is formed from A as in FADDD. 
FLOATING SUBTRACT MAGNITUDE 


mnemonic address code no. of words 


FSBM Xx 10 2 


The absolute magnitude of the floating point number in X and 


X + 1 is subtracted from FAC. The result replaces FAC. 
DEFERRED FLOATING SUBTRACT MAGNITUDE 


mnemonic address code no. of words 


FS BMD A 110 2 


Deferred subtract magnitude. Effective address is formed from 


A as in FADDD. 
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2.2.3 TYPE III COMMANDS - UNARY COMMANDS 


FLOATING ABSOLUTE VALUE 


mnemonic address code ' no. of words 
FABS none 14 1 


The absolute value of the FAC replaces the FAC, i.e. | FAC | 


replaces FAC. 
FLOATING SQUARE 


mnemonic address code no. of words 


FASQ none 15 1 


The square of FAC is returned in FAC. This instruction requires 
that the mantissa of FAC be normalized prior to execution as in 


type II instructions (see 2.2.2). 
FLOATING NORMALIZE 


mnemonic address code no. of words 


FNOR none 16 1 


The contents of FAC are normalized and replace FAC. This in- 
struction can cause exponent overflow or underflow in which case FAC 
will contain the largest possible negative floating point number or 


all zeros respectively and FXFLG will be set non-zero. A successful 


normalize will set FXFLG to zero. 


FLOATING NEGATIVE VALUE 


mnemonic address code no. of words 


FNEG none 17 1 


The contents of FACHI and FACLO are twos complemented, i.e. -FAC 


replaces FAC. 
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2.2.4 TYPE IV COMMANDS - INDEX COMMANDS 
LOAD INDEX 
mnemonic address code no. of words 
FLDX I 27 2 


The pseudo-index is loaded with the 16 bit contents of location I. 


STORE INDEX 


mnemonic address code no. of words 
FSTX I 30 2 


The 16 bit pseudo-index is stored into location I. 


2.2.5 TYPE V COMMANDS - CONDITIONALS 


These commands allow the program to alter the path of control which the 
interpreter is following based on the results of certain tests. The lo- 
cation to which the interpreter is caused to transfer must contain a valid 
floating point command. If the interpreter should encounter an invalid com- 
mand at any time during execution, it will come to a halt with the address 
of the illegal command displayed in the MB register on the front panel. This 


is the only halt in the program. 
JUMP UNCONDITIONAL 


mnemonic address code no. of words 


FJMP C 20 . 2 


Unconditional jump. The interpreter will take the next command 


from location C and continue from there. 
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JUMP IF AC POSITIVE 
mnemonic address code no. of words 
FJAP C 21 2 


If FAC is positive or zero, the interpreter takes the next 
command from location C. Otherwise, the interpreter continues with 


the command following the FJAP command. 


JUMP IF AC ZERO 


mnemonic address code no. of words 


FJAZ C 22 2 


If FAC is 0, the interpreter will take the next command from 
location C. Otherwise, the interpreter continues with the command 
following the FJAZ command. Note; The interpreter tests only FACHI 
for zero. FAC may be non-zero and FACHI = 0 only if the number in 
FAC is not normalized. This condition cannot be created by the inter- 


nnormalized numbers into his 


calculations (see 2.2.2). 
JUMP IF AC NEGATIVE 


mnemonic address code no. of words 


FJAN C 23 2 


If FAC is negative, the interpreter will take the next command 
from location C. Otherwise, the interpreter continues with the com- 


mand following the FJAN command. 
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JUMP IF EXPONENT OVERFLOW (OR UNDERFLOW) 


mnemonic address code no. of words 


FJEV C 24 zZ 


If FXFLG is non-zero, the interpreter will take the next command 

from location C and set FXFLG to zero. Otherwise, the interpreter will 
continue with the command following the FJEV command. The FJEV command 
is used to detect the occurrence of either exponent overflow or exponent 
underflow resulting from the execution of the last preceding Type II command 
or FSTA, FNOR, or FASQ. If desired, the type of overflow may be detected 
by an FJAZ command at location C, since exponent underflow returns FAC=0, 
and exponent overflow returns the largest number (+ or -) in FAC. 
JUMP IF DIVIDE CHECK 

mnemonic address code no. of words 

FJDC C 25 2 
If FDFLG is non-zero, the interpreter will take the next command 

from location C and set FDFLG to zero. Otherwise, the interpreter con- 
tinues with the command following the FJDC command. The FJDC command is 
used to detect the occurrence of divide check during execution of the last 
previous FDIV or FDIVD command. If desired, one may test whether the con- 
dition occurred because the divisor was 0 or not normalized by checking the 
divisor with an FLDA and FJAZ instruction at location C. 
JUMP IF INCREMENTED INDEX NOT ZERO 

mnemonic address code no. of words 


FJIX C 26 : 2 


The pseudo-index (FINDX) is incremented by one, and if the result is 


non-zero, the interpreter takes the next command from location C. 
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If the result is' 0, the interpreter continues with the command follow- 
ing the FJIX instruction. The pseudo-index will contain the incremented 


value whether or not the jump occurs. 
2.2.6 TYPE VI COMMAND - EXIT 
EXIT FLOATING INTERPRETER 


mnemonic address code no. of words 


FEXT none 0 1 


This command causes the interpreter to return control to the 
user at the location immediately following the FEXT. None of the 
internal registers or flags are altered by either the FEXT or entering 


the package. The AO is returned in the ADD state. 
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CHAPTER THREE 


DATA CONVERS ION 


Introduction: 


Two conversion routines are provided; one to convert from floating point 
to character, the other to convert from character to floating point. Both 
conversion routines are core to core operations rather than being bound to a 
particular 1/0 device (that is, characters are fetched from and stored into 
memory). For added flexibility, all characters are referenced with an index 
into a character set table called @FCST, which initially contains 8-bit 
ASCII codes. Changing the character set for a specific I/0 device can easily 


be accomplished by changing the character codes in @FCST. 


External floating point format is expressed as a mantissa or fraction 
portion and a power of ten by which the mantissa is multiplied. This is 


written as tn.nnnnnn E +nn, where n is a decimal digit. The number to the 
Lm ; ——— 


mantissa exponent 
right of the E is the power of ten by which the mantissa is multiplied. Thus, 
-~3.527614E+03 is -3.527614*103 or -3527.614. The floating point number 


+172.100123E-02 is +172.100123*1072 or +1.72100123. 
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3.2 Floating Point to Character Conversion: 


NAME; 


SUBROUTINES CALLED: 


ALTERED REGISTERS & FLAGS: 


CALLING SEQUENCE: 


ARGUMENTS : 


FUNCTION: 


@SFC 


SSFI, @FXC 
FAC, FIM, FXFLG, FDFLG 


JU @SFC 
WRD el-1 
WRD e2 
WRD e3 


return 


el is the address of the location into 
which the first output character is 


to be stored. 


e2 is the address of the two word float- 
ing point argument. The argument need 
not be normalized but the magnitude must 
be zero or in the range (27129 , +128) 
(in decimal this is 1.469367E-39 to 
1.701411E+38) 


e3 is the address of the error return. 


Converts a signed two word floating 
point argument to a string of thirteen 
characters, stored one character per 
word, right justified starting in lo- 
cation el. The character string is of 
the format 


“+ 
i 


’ n.nnnnnn z|+/on 
= 

Ln 

where n is character representation of a 


decimal digit. 


ERRORS : 


NOTE: 


LENGTH: 


Description of Algorithm: 
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If normalization of the floating point 
argument caused either exponent over- 
flow or underflow, an * is stored rather 
than a leading + or - sign, and when 
conversion is completed, control returns 
to e3. The * can be considered a - sign. 
An argument resulting in overflow con- 
verts to * 1.701411 E + 38. An argu- 
ment resulting in underflow converts to 
*0.000000 E + 00. 


The magnitude of the three smallest nor- 
malized non-zero floating point numbers 
are converted to one of the character 
strings +1.469367E-39 or +1.469368E-39. 
These two character strings cannot be 
converted back to a floating point number. 
The smallest character string which can 
successfully be converted to a floating 
point number is +1.469369E-39. Therefore, 
if the user converts any one of these 
three numbers to a string of characters, 
he should be aware that he cannot suc- 
cessfully convert the string back to 


a‘ floating point number. 


306g (198,9) locations 


The sign of the floating point argument is stored, the argument is then 


normalized, and the absolute value is taken and used for conversion. If nor- 


malization caused either exponent overflow or underflow, the error return is 


taken when conversion is completed and an asterisk (which may be considered 


asa' 


-'), is stored rather than a leading sign. 
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ee, 
Since the output character string is of the form /+ n.nnnnnnEtnn, the 


floating point argument is first manipulated to make ae than or 
equal to one and less than ten. (If the floating point argument is exactly 
zero, this portion of the algorithm is bypassed.) Making 1£ floating 

point argument <10 is accomplished by first checking if it is>1. If it 
is not, it is multiplied by the largest possible power of ten (1038) and, 
if necessary, it is multiplied once again by ten to force it21. The ar- 
gument is then checked for< 10. If it is not<10, it is forced so by 
dividing by the largest power of ten, which is less than the argument. The 
powers of ten used in multiplying and dividing the argument to force its 


value to be between one and ten are used to form the exponent portion of 


the character string. 


With the floating point argument (stored in FAC) now21 and <10, the 
mantissa portion of the character string can be formed. FACHI, FACLO is 
treated as a double precision mixed number with FACXP showing the position 
of the binary point. FACHI, FACLO is left shifted (with overflow bits shifted 
into a 3rd word) until the binary point immediately precedes bit 15 of FACHI. 
The overflow word is then converted to character and stored as the first digit 
of the mantissa, immediately followed by a decimal point. The fraction por- 
tion of the mantissa is formed by successively multiplying FACHI, FACLO by 
1059 and storing the most significant word of the 3 word product. The ex- 


ponent is then converted and stored, preceded by an E and either a + or - sign. 
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3.3 Character to Floating Point: 
NAME: @SCF 
SUBROUTINES CALLED: SSFI, @GFXC 
ALTERED REGISTERS & FLAGS: FAC, FIM, FXFLG, FDFLG 
CALLING SEQUENCE: JU @SCF 
WRD el-1 
WRD e2 
WRD e3 
return 
ARGUMENTS : @l is the address of the first character 


in the string to be converted. The char- 
acter string should be stored one charac- 
ter per word right justified in the format 


\ gi 
[+]( nfn}]i [4...] Cd] 


n [news] s [oie] 


Le Tis) wee } E 


The notational conventions are: 


1. nis a decimal digit 

2. A&A is a space 

3. g is a delimiter 

4. braces [] contain optional items 
which may or may not be included. 

5. braekets {{ contain alternate items 
where one and only one of the items 
must.be included. 

6. ellipses ... denote permissible 
repetition of the preceding iten. 


The string is treated as follows: 


1. If there is no sign, it is treated as +. 
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2. If the leading sign is * or -, it 
is treated as -. 

3. If there is no decimal point, it is 
assumed to follow the last mantissa 
digit. 

4. Characters are processed up to and 
including the first 2, or 1349 char- 


acters have been processed. 


e2 is the address where the two word 


floating point answer is stored. 
e3 is the address of the error return. 


FUNCTION: Converts a string of decimal characters 
to a two word normalized floating point 
answer. The two word normalized float- 
ing point answer is returned in registers 
AX (MSH), AY (LSH), and is stored in lo- 
cation e2 (MSH) and e2+1 (LSH). The AO 


is returned in the ADD state. 


ERRORS: A scan error occurs if the character 
string is illegally formed. Location 
@SCF+3 is set to zero and control im- 


mediately returns to e3. 


An overflow error occurs if the charac- 

ter string contains more than 1049 man- 
tissa digits (discounting leading zeros) 

or if the magnitude of the number is out~- 
Side the range 1.469369E-39 to 1.701411E+38. 
Location @SCF+3 is set to one and control 


immediately returns to e3. 


Whenever control returns to e3, the AO 


is in the ADD state. 
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NOTES: See NOTES under @SFC. 
LENGTH: 406. (26219) locations 


Description of Algorithm 


The mantissa portion of the character string is converted to a 
double precision integer by multiplying the answer by 1049 and adding 
in the latest digit. This double precision mantissa is then converted 
to a normalized floating point number. A count of the number of digits 
to the right of the decimal point is kept and, after the exponent portion 
of the character string has been converted, this digit count is subtracted 
from it to obtain the final exponent. The magnitude of the final exponent 
is used as an index into the positive floating point powers of ten table 
(see 3.6). The floating point number obtained from the mantissa portion 
of the character string is then multiplied (if the final exponent was 
positive) or divided (if the final exponent was negative) by this power 
of ten to form the final floating point answer. If there was a leading 
minus sign or asterisk, the floating point answer is two's complemented 


before return. 


3.4 Common Tables & Routines; 


The conversion routines @SFC and @SCF reference a common routine called 
@FxC, which has four entry points. @FXC occupies a total of 213. (13949) 
locations. Since @FXC is common to both @SFC and @SCF, it need appear only 
once if the conversion routines are used together. In the discussion of @FXC 
which follows, each of the four entry points is treated separately for the 


sake of clarity. 
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3.5 Character Set Table: 
NAME: @FCST 
FUNCTION: Common external character set table for 


floating point data conversion routines. 


The table is ordered as follows: 


Location Contents 

@FCST code for zero 

@FCST+1 code for nine 

@FCST+2 code for + 

@FCST+3 code for - 

@FCST+4 code for * 

@FCST+5 code for . 

@FCST+6 code for E 

@FCST+7 code for space 

@FCST+10 code for delimiter 
NOTES: The standard table is in full 8-bit ASCII. 


The delimiter character at @FCST+102 is a 


carriage return and may be changed if 


desired. 


The entire table may be replaced with a 
different character set provided that the 
numeric codes in the new set are sequential 
and the code for zero (0) is less than the 
code for nine. No code may occupy more 
than 15 bits. 


LENGTH: 129 (1010) locations 
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3.6 Floating Point Powers of Ten Table 
NAME: @GFPT 
FUNCTION: Common floating point positive powers of 


10jq table for floating point data con- 
version routines. Each floating point 
power occupies two locations in the table. 


The table is organized as follows: 


@FPT: WRD 45473,46777 constant for 1038 
WRD 74136,160773 constant for 1037 
WRD 40000, 201 constant for 100 
NOTE: | @FPT is located at @FCST+12, 


LENGTH: 116g (78,9) locations 


3.7 Left Shift FAC: 


NAME: 
SUBROUTINES CALLED: 
ALTERED REGISTERS & FLAGS: 


CALLING SEQUENCE: 


ARGUMENTS : 


FUNCTION: 


ERRORS : 
NOTES : 


LENGTH: 
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@LSHF 
none 
none 


Load AX with the negative shift count 
JU  @LSHF | 


return 


Register AX contains minus the number of 


places to left shift FACHI, FACLO 


Performs double precision left shift of 
FACHI, FACLO. On return, the shifted 
result is in AX (MSH), AY (LSH). Any 
carry out of MSH is found in location 


@DIG+1. 
none detected 
@LSHF=@FCST + 3068 


@DIG+1=G@LSHF + 14g 
23 (19) 6) locations 


3.8 


Multiply FAC by Ten: 


NAME: 
SUBROUTINES CALLED: 


ALTERED REGISTERS & FLAGS; 
CALLING SEQUENCE: 
ARGUMENTS ; 


FUNCTION: 


ERRORS : 


NOTES; 


LENGTH: 
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@10X 
@LSHF 


FAC 
JU @10x 
n/a 


Performs unsigned multiplication of 
FACHI, FACLO by 10,9. The most sig- 
nificant word of the three word product 
is returned in AY. The second and 
third words of the product are found 
in FACHI, FACLO respectively. 


n/a 
@10X=@FCST+650 


40, (3216 ) locations 
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CHAPTER FOUR 


EXTENDED COMMANDS 


Introduction: 


In addition to the basic floating point interpreter, a set of mathematical 
functions is supplied which can be invoked by a command in the same line with 
the basic commands. These functions also call the floating interpreter and 
since the interpreter has already been entered at this point, a push-down 
scheme is supplied to allow recursive calls such as this. The push-down list 


will accomodate recursive calls up to seven levels. 


It should be noted that all pseudo registers - the floating accumulator, 


common to all levels of the recursion. In other words, if an extended func- 
tion which calls the interpreter recursively is invoked by a command, these 
registers and/or flags may be altered. Information detailing such factors 


is supplied in the documentation accompanying the individual package. 


The push-down scheme and command code structure is tailored so that the 
user may easily add his own functions. The procedure for doing this is 


described in section 5.3. 


The mathematical subroutines which are supplied with the extended package 
are SINE, COSINE, ARC TANGENT, LOG,, EXPONENTIAL, and SQUARE ROOT. The com- 
mands associated with these are FSIN, FCOS, FAIN, FLNE, FEXP, and FSQT (codes 
31, 32, 33, 34, 35, and 36) respectively. They each perform the desired 
function on the contents of the floating pseudo-accumulator and return the 
results in the same register. Errors which can result, such as attempting 


to take the square root or log of a negative number, are flagged by the 
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routines in internal locations not accessible in interpretive mode, i.e. 
cannot be tested with an interpreter command. An error trap routine is 
available which will handle these and other errors when they occur (see 


Appendix E). 


In the writeups that follow, FAC is the floating pseudo-accumulator, 
FIM is the temporary floating pseudo-accumulator, FDFLG is the divide check 


flag, FXFLG is the exponent overflow flag, and FINDX is the pseudo-index. 


4.2 Sine, Cosine: 


COMMAND: 


FUNCTION: 


ERRORS ; 
ALTERED REGISTERS & FLAGS: 


METHOD: 
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FSIN (code 31), FCOS (code 32) 


a. FSIN - calculates the SINE of the 
contents of FAC which is assumed to be 
a radian argument and replaces FAC 


with the result. 


b. FCOS - calculates the COSINE of 
the contents of FAC which is assumed to 
be a radian argument and replaces FAC 


with the result. 
none 
FAC, FIM, FXFLG 


For FCOS, the absolute value of FAC is 
subtracted from **/2 (=1.570796) and the 
SINE of the result is taken. 


For FSIN, the argument (FAC) is first 
multiplied by 2/]] to convert it into 
units of a quarter circle, and the result 
is checked for its absolute magnitude 
being less than one. If so, it is a 
first quadrant quantity and the procedure 
continues with the series calculation 
described later. If the magnitude of the 
result is greater than or equal to one, 
its sign is saved, it is forced positive, 
and the integer portion is shifted out - 
leaving a positive fraction (referred to 
as Y in the following). The last two 
bits of the integer portion and the sign 
are used to determine which quadrant the 


original argument was in and the quantity 


Y is altered as follows: 
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sign last two bits 4 quadrant 

+ 00 Y ®Y I 

+ 01 1-Y »Y II 

+ 10 -Y?Y III 

+ dei alr Yay IV 

- 00 =Y =X IV 

- 01 LP YY : IIL 

- 10 YoY II 

= 11 Lay Sy I 


This new value of Y is then treated as a fraction and is normalized. 


The series used to calculate the sine is basically a 5 term Chebyshev 
economized polynomial approximation of a 6 term McLaurin series for sin | 2 r 
The coefficients are further "adapted" to allow the series to be calculated 
with one less multiplication than would be the case for a standard polynomial 


evaluation procedure. This results in the sine being calculated as follows: 
sin = ((Z -Yt Ay) * Z + Az) * A, * Y 
where 
Z= (Y+Aj) *Y +A) 
and 
Ag = -14.93104811 
Ay = 739.74079011 
Ay = +367.8139482 
A3 = +23410.00773 


A, = +0.0001514440767 


Accuracy is 6 + significant decimal digits for arguments in the first 
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quadrant () rac(Z%. Accuracy loss is about two thirds of a decimal 


digit for each complete rotation, i.e. if Tn <[racl 4277 (+1), the ac- 


curacy is about 6 - en decimal digits. 


4.3 


Arc Tangent: 


COMMAND : 


FUNCTION: 


ERRORS ; 


ALTERED REGISTERS & FLAGS: 


METHOD: 


4-6 
FLOATING POINT MANUAL 
71-44-001 


FATN (code 33) 


The arc tangent of the contents of FAC 
replace FAC. The result is in radians 


and lies in the range ( Bel +O, 
none 


FAC, FIM, FDFLG, FXFLG 


The argument (FAC) is checked for its 

absolute magnitude being greater than or 
equal to one. If so, a flag is set and 
the reciprocal of the argument is taken 


and replaces FAC. 


The arc tangent of the quantity in FAC 


is then approximated by 


(An + Ay X24 A5X4) 
Z = ATAN X= X* (Bo + BLX2+ BoX4) 


where X is the argument and 


Ay = 0.6402481953 
A, = 0.4229908144 
Ay = 0.0264694361 
Bg = 0.6402487022 
By = 0.6363779373 
By = 0.1108328778 


If the flag was set by the initial check, 
the value Z is checked for + or -. If Z 
is +,(+ - Z)replaces Z. If Z is -, 

(- +z) replaces Z. (This is effected 
by subtracting Z from + or - T depending 
on the sign of Z.) 

If the flag was not set by the initial 
check, the value Z is not altered. Accur- 


acy is 6+ significant decimal digits for all 
arguments. 


4.4 


Natural Log: 


COMMAND: 


FUNCTION: 


ERRORS : 


ALTERED REGISTERS & FLAGS; 


METHOD: 
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FLNE (code 34) 


The natural log of the contents of FAC 
replace FAC. 


If FAC is negative, a flag (FNLNF) is set, 
FAC is forced positive, and the natural 


log taken. 


FAC, FIM, FXFLG, FNLNF (FPLNE+4) 


The quantity in FAC is 
Z=xX- 2! where .54x<1 and I is an 


integer. 


In Z = 1n [x°° 24] 
In X + I 1n2 


The quantity In X is approximated by the 
polynomial. 


In X= In A-Z (¥ + Y3/3 + Y?/5 + ¥//7) 


which is a Taylor series evaluated at A 


The product [I ln 2] is added to ln X, and 
the sum is left in FAC. 


A = 0.70710678 
In A = 0.34657359 
In 2 = 0.69314718 


Accuracy is 6+ significant decimal digits 
except for .9044241.110. In the latter 


range, accuracy decreases as Z~1l. 
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4.5 Exponential: 
COMMAND: FEXP (code 35) 
FUNCTION: The exponential of the contents of FAC 


replace FAC. (FAC = eFAC) 


ERRORS : If the result is going to be out of 
range, i.e. if FAC 88.722, a flag 
(FEXOF) is set. If FAC was negative, 
zero is left in FAC. If it was positive, 


the largest positive number is left. 


ALTERED REGISTERS & FLAGS: FAC, FIM, FDFLG, FXFLG, FEXOF (FPEXP+1) 


METHOD: eX = 2% log e 
where I is the integer portion 
and F is the fractional portion of X log, e 


Multiplication by 2! is computed by the 


continued fraction: 


A 
BtFt+t+C 
F+D 
F 
where 


A = -34.624680982 
B = -17.312340491 
© = 104.0684491 
D = 20.813689813 
logy e = 1.442695041 


Accuracy is 6+ significant decimal digits 
for {x|g10. Accuracy decreases slowly as 
|x| becomes large until at {x[=ss, the ac- 
curacy is 5+ significant decimal digits. 


4.6 


Square Root: 


COMMAND: 


FUNCTION: 


ERRORS: 


ALTERED REGISTERS & FLAGS; 


METHOD: 
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FSQT (code 36) 
The square root of FAC! replaces FAC. 


If FAC is negative, it is forced posi- 
tive, and FSFLG (internal to the square 
root routine) is set non-zero. If FAC 


is positive, FSFLG is set to zero. 


FAC, FSFLG (=FPSQT + 6) 


After FAC is forced positive and FSFLG 
is determined, the exponent of the re- 
sult is determined by dividing FACXP 

by two (by shifting right once) and 

adding 100g to preserve the excess 200g 
notation. If the original exponent was 
odd, the shifted FACXP is increased by 


one; otherwise, it is left alone. If 


and FACLO are shifted left once. Since 
the algorithm treats FACHI and FACLO as 
a 32 bit positive fraction with the 
binary point to the left of bit 15 of 
FACHI, the fact that the left shift will 
set the sign bit (bit 15) of FACHI does 


not matter. 


The algorithm then proceeds to determine 
a fourteen bit first approximation to 
the square root by a method based:»on the 
fact that N* is the sum of the first N 
odd numbers. This method also leaves as 
a "remainder" the difference between the 


square of the approximation and the original 
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number. This remainder and the initial 
approximation are then used for one 
Newton-Raphson iteration which completes 
the square root using the single pre- 
cision divide entry (FSDVD) of the 
floating point package. 


Accuracy is 6+ significant decimal 


digits for all input arguments. 
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CHAPTER FIVE 


NON-INTERPRETIVE MODE USAGE 


Introduction: 


Certain sections of the floating point interpreter are directly 
accessible to the user without the need to supply commands. These sec- 
tions may be invoked by a JU SUBR instruction and, after the operation 
is completed, will return control to the instruction following the jump. 
In order to use these routines successfully, it is necessary to know that 
in addition to the pseudo-accumulator (FACHI, FACLO, and FACXP) there is 
a "temporary" accumulator (FTMHI, FIMLO, and FTMXP) which is used to 
contain the floating argument of a Type II command during the execution 
of oe operation (see 1.3). This temporary pseudo-accumulator, referred 
to as FIM, is loaded in the same manner as FAC (see FLDA instruction in 
2.2.1). If the user desires to access the routines described in this sec- 
tion, he may need to load FIM in addition to FAC for those routines that 


operate on both accumulators. 


These sections wili be described as subroutines since they are essen- 
tially used in this manner when accessed directly. When the floating in- 


terpreter resides in memory, all of these subroutines also lie in memory. 


A "command equate" source tape is included in the software package. 
This tape defines these subroutines as well as the various pseudo-registers 


and locations associated with them (see chapter 6). 
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5.2 Subroutines: 
5.2.1 Double Precision Fixed Point Add 
NAME: FDAD 
CALLING SEQUENCE; JU FDAD 
INPUT: FACHI, FACLO; FIMHI, FIMLO; AO must be in 
ADD state. 
FUNCTION: FACHI, FACLO and FIMHI, FIMLO are treated 


as signed double precision numbers and 
added. The result of the addition ap- 
pears in FACHI, FACLO. FTMHI, FIMLO are 
left unchanged. 


If arithmetic overflow occurred (two 
numbers of like sign are added and the 
result has opposite sign), the link 
will be set to 1. If no arithmetic 
overflow occurred, the link will be 


zero. 


The AO is in the ADD state upon return. 


NOTES ; It is possible to generate the maximum 
negative number (FACHI = 100000, FACLO = 
000000,), which is not considered a case 
of arithmetic overflow; and so the link 


will not be set. 
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5.2.2 Double Precision Fixed Point Multiply 


NAME: 
CALLING SEQUENCE: 


INPUT: 


FUNCTION: 


NOTES ; 


FDMPY 
JU  FDMPY 


FACHI, FACLO; FTMHI, FIMLO 

AX must be set to the value in FTMHI 
AY must be set to the value in FTMLO 
The AO must be in the ADD state 


FACHI, FACLO and FTMHI, FIMLO are 
treated as signed double precision 
numbers and are multiplied. The high- 
order 30 bits of the 62 bit product 
are returned, right justified, in 
FACHI, FACLO. The value in FTMHI, 
FTMLO is unchanged. 


The AO is in the ADD state upon return. 


The 30 bit product is inaccurate in the 
right-most two bits. If FACHI, FACLO 
and FTMHI, FIMLO are each considered as 
a double precision fraction with its 
binary point immediately to the right 
of the sign, i.e. between bits 14 and 
15 of the high-order word, the binary 
point of the product will be shifted 
right once so that it is between bits 
13 and 14 of FACHI. 
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5.2.3 Double Precision Fixed Point Divide 


NAME: | FDDIV 
CALLING SEQUENCE: JU FDDIV 
INPUT: FACHI, FACLO; FTMHI, FTMLO 


AX must be set to the value in FIMHI 
AY must be set to the value in FTMLO 
The AO must be in the ADD state 


FUNCTION: FACHI, FACLO and FIMHI, FIMLO are 
treated as signed double precision 
numbers, and the former is divided 
by the latter. The quotient appears 
in FACHI, FACLO. The value in FTMHI, 
FTMLO has been destroyed. 


The quotient will be 30 bits in FACHI, 
FACLO with the binary point displayed 
one position to the right in the same 


way as explained in the note for FDMPY. 


The absolute magnitude of FTMHI, FTMLO 
must have bit 14 of FIMHI set for the 
divide to take place. If this con- 

dition is not satisfied, divide check 


will occur. 


The AO is in the ADD state upon return. 


NOTES: ' The rightmost three bits of the quotient 
are inaccurate. Divide check causes 
FACHI, FACLO to be set to a large double 
precision number of the sign which would 
result if the divide could take place 
(FACHI, FACLO = 077777, 177400 or 100000, 
000400 for + and - respectively); also, 
FDFLG is set non-zero. A successful 


divide sets FDFLG to zero. 
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5.2.4 Single Precision Divide 
NAME: | FSDVD 
CALLING SEQUENCE: JU FSDVD : 
INPUT: , . AX = high order dividend ; must be a posi- 


4 
FLODV = low order dividend ‘?V % Pit 
‘double preci- 
sion number (see 
below) 


AY = negative divisor 


The AO must be in the ADD state. 


FUNCTION: This is an inner loop which, if used 
correctly, can be invoked to supply an 
unsigned single precision divide. The 
quotient is incomplete in the sense that 
it is right shifted and truncated upon 


return. 


To obtain a complete single precision 
unsigned divide, the following procedure 
may be used. First, load AX and the lo- 
cation FLODV with a valid two word posi- 
tive product (bits 14 and 15 of AX must 
be zero). Then load AY with the positive 
single precision divisor and twos comple- 
ment it. The following code will then 
perform the divide: 


JU FSDVD 3 LNCOMPLETE QUOTIENT IN TRP 
RRC AO,L1,0 3;GET LAST BIT OF QUOTIENT 
This code may be SFM NOT LNK 3; UPDATE 
eliminated if the 
remainder is to be sca n0see eee 
disregarded. NOP 3IN AX 


RR - TRP,L1,AY 3TRUE QUOTIENT IN AY 
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Note that the incomplete quotient is in 
the TRP register on return from FSDVD. 
The AO is in the ADD state upon return. 


If either the link is set or AY (the 
final quotient) is negative following 
this code, divide check has occurred. 
This means that the high-order portion 
of twice the dividend was greater than 
or equal to the divisor, and the quotient 


is incorrect. 


NOTES : No flag is set if divide check occurs. 


5.2.5 Floating Point Normalize 


NAME: 
CALLING SEQUENCE: 
INPUT: 


FUNCTION: 


1) 


or, 
2) 
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FNORM 
JU FNORM 
FACHI, FACLO, FACXP 


Same as FNOR command (see 2.2.3), in- 
cluding the setting of FAC and FXFLG 
should exponent overflow or underflow 


occur. 


The advantage of the accessibility of 
this routine lies mainly in the saving 
of time. For instance, to convert a 
single precision integer value to 
floating point, the following two 
methods could be used. (Assume the 
integer is in AX, and the floating 


equivalent is wanted in location X.) 


RM AX, FACHI 
ZM FACLO 
MRI 217, AX 
RM AX, FACXP 
JU FNORM 

JU SSFI 
FSTA X 

FEXT 

RM AX, FACHI 
2M FACLO 
MRI 217, AX 
RM AX, FACXP 
JU SSFI 

FNOR 

FSTA X 
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Version 1) takes one more location in 


core and saves about 80 machine cycles. 


NOTE: The AO may not be in the ADD state upon 


return. 


5.2.6 Negation and Store 
NAME: 

CALLING SEQUENCE: 

INPUT: 


FUNCTION: 


NAME: 


CALLING SEQUENCE: 


FUNCTION: 


NAME: 
CALLING SEQUENCE: 
INPUT: 


FUNCTION: 


NAME: 
CALLING SEQUENCE: 
INPUT: 


FUNCTION: 
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FACMP, FACMA 
JU FACMP or JU FACMA 
FACHI, FACLO or AX, AY 


a) FACMP - replaces FACHI, FACLO with 
its two's complement. Result is 
also returned in AX, AY. 

b) FACMA - replaces FACHI, FACLO with 
the two's complement of the double 
precision number in AX, AY. Result 
is also returned in AX, AY. 


FICMP, FICMA 


JU FICMP or JU FTCMA 


a) FTCMP - replaces FIMHI, FTMLO with 
its two's complement. Result is 
also returned in AX, AY. 

b) FTCMA -— replaces FTMHI, FIMLO with 
the two's complement of the double 
precision number in AX, AY. Result 


is also returned in AX, AY. 
FASAX 
JU FASAX 
AX, AY 


Stores AX into. FACHI and AY into FACLO 


FTSAX 
JU FTSAX 
AX, AY 


Stores AX into FTMHI and AY into FTIMLO 


~10 : 
oar FLOATING POINT MANUAL 


71-44-001 
5.2.7 Generate Zero or Largest Number 
NAME: FOF AC 
CALLING SEQUENCE: JU FOFAC 
INPUT: none 
FUNCTION: sets FACHI, FACLO and FACXP to zero 
also returns AX and AY = 0 
NAME: FCMAX 
CALLING SEQUENCE: JU FCMAX 
INPUT: | FACHI 
FUNCTION: FACHI, FACLO, FACXP will be set to the 


maximum possible floating point number 


of the original sign of FACHI. 


1) If FACHI< 0, this routine sets 


FACHI = 100000 
FACLO = 000400 
FACXP = 000377 


2) If FACHI=0, this routine sets 


FACHI = 077777 
FACLO = 177400 
FACXP = 000377 


Upon return, AX and AY will be equal to 
the value stored in FACHI and FACLO re- 


spectively. 
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5.2.8 Floating Arithmetic Right Shift 
NAME: FARSN 
CALLING SEQUENCE: JU FARSN 
INPUT: AX, AY, FARSC 
FUNCTION: | This routine arithmetically right shifts 


the double precision number in AX, AY 
by the number of places indicated by 
-~FARSC e 


NOTES ; FARSC must be set to a negative count 
before calling FARSN. 


oF 
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5.2.9 Other Notes on Non-Interpretive Usage 


1) The pseudo-index is kept in location FINDX and may be set by the 
user without using an FLDX command (see 2.2.4) by simply storing 
the desired value via a RM R,FINDX where R is a register con- 


taining the index value. This, as with FNORM, is a time saver. 


2) The two flags, FXFLG and FDFLG, are in locations defined by their 
names, and can be checked (or cleared) in non-interpretive mode 


to save time. 


3) The usage of the locations FTBLE, FARGD, FETCH, and FMASK which 


are on the command equate source tape is described in 5.3. 


User Generated Extended Functions; 


If the user desires to add functions of his own to the extended package, 


the procedure is quite easy as outlined below. 


The extended package as delivered uses command codes 00-36, inclusive 
and 101, through 110. inclusive. There are available codes of 379-77, inclu- 


sive which the user may assign to his own functions. 


User functions may be of two types - invoked by one word commands or 
invoked by two word commands where the second word is an argument address 
or value. If deferred mode addressing ie Gesared as ba option for the same 
function, it must be accomplished by user code. Setting bit 6 of the command 
code to attempt deferred addressing will cause the floating interpreter to 


take the error halt. 


Suppose the command name used to invoke the function is to be FFCN 
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assigned to code 379. 


Step 1) 


Step 2) 


Step 3) 


tep 4) 


Using the Source Text Editor, add the statement FFCN = 37 


to the Command Equate Tape (see operating instructions). 


The user code which accomplishes the function must have 


the following code just before the END statement 


LOC FTBLE + FFCN 
WRD ENTRY 
where ENTRY is the location at which the user function 


begins execution. 


The last instruction executed by the user function must re- 
turn control to FGET, usually via a JU FGET. Remember that 
when the user function is invoked by a command, the interpreter 


passes control to the user function. The JU FGET returns control 
to the interpreter. 

If the function the user is generating needs the fioating point 
capability supplied by the interpreter, the user function may 


call the interpreter followed by a list of commands to ac- 


complish the task subject to the following restrictions: 


a) The command name corresponding to the function it- 


self (in this case, FFCN) may not be used. 


b) Commands which cause the interpreter to be called re- 
cursively may be used so long as care is taken not to 
exceed seven levels of recursion in total (see 4.1) 
(remember that the function being coded is at least 
at level 1 during its execution, and if it calls the 
interpreter, all commands in the list are at least at 
level 2). 


c) No function invoked by a command may have in its code 
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a call to the interpreter whose command string con- 
tains the command name corresponding to the function 
itself. This is an indirect violation of restriction 


a) above. 


Step 5) Assemble the function using the Command Equate Tape for 


pass 1 as explained under operating instructions. 


Notes: If the function being generated is invoked by a two word com- 
mand whose second word is an argument, one and only one of the follow- 


ing steps must occur during its execution. 


a) JU FARGD 
This fetches the contents of the location following the command into 
register AX. 
b) JU FETCH 
This calls FARGD and uses the contents of the location following the 
command as an address to fetch a floating point argument which is 
placed in FTM. Also, AX and AY will be set to the value in FTMHI 
and FTMLO respectively upon return. 
c) ZM FMASK 

JU FETCH 
This causes deferred fetching of a floating point argument. The con- 
tents of the location following the command is used as an address of 
another location which is incremented twice to form the addresses of 


the floating argument which is loaded into FIM and AX, AY as in b). 


Examples: 
1) FCSX is to be the command name, 1 word, code 37. When in- 


voked it is to take the COSINE of the SINE of the value in FAC. Assume 


Step 1 has been accomplished by adding the statement FCSX=37 to the command 
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equate tape. This function may be accomplished by the following code: 


FCS: JU SSFI . senter floating interpreter 
FSIN 3sin of FAC | 
FCOS ;cos of FAC 
FEXT ;exit interpreter 
JU FGET ;return to interpreter 
LOC FIBLE + FCSX 
WRD FCS 


END 
When this is assembled with the Command Equate Tape resulting 
from step 1 and loaded with SSFI and the SINE, COSINE routine, the 


user may now call the routine in the floating interpretive mode as 


follows: 


FCSX 


FEXT 
Note: The routine must be assembled with the new command equate tape. 
2) FMCS is to be the command name, 2 words, assigned to code 40,. 
When invoked, it is to take the SIN of the COSINE of the value in FAC 
and set the sign of the result to the sign of the floating point ar- 
gument whose address is the second word of the command. Step 1 re- 
quires the new command to be added to the command equate tape. This 


function could be coded as follows: 


3-16 
FLOATING POINT MANUAL 


71-44-001 

FMC: JU FETCH 3fetch arg to AX,AY 

RMI AX, 0 3;save MSH arg (sign of arg) 

JU SSFI s;enter floating interpreter 

FCOS 3cos of FAC 

FSIN s3sin FAC 

FABS jabs value of FAC 

FEXT ;exit floating interpreter 

MR FMC + 3,AX ;get sign of arg 

JC AX, GEZ, FGET j;plus, exit 

JU FACMP sminus, comp FAC 

JU FGET ;return to interpreter 


LOC FIBLE + FMCS 
WRD FMC 


END 


When assembled with the Command Equate Tape and loaded with SSFI, 
and the SINE, COSINE routine, it may be invoked by another routine 


assembled with the Command Equate Tape via 


JU SSFI 
FMCS X 
FEXT 


3) FMCSD is to be the command which does the same thing as 
FMCS, only using deferred mode addressing for the argument. FMCSD 
must be assigned a different code - say 41... Both FMCSD and FMCS 


may be coded in the same routine as follows, assuming their equates 


have been added to the Command Equate Tape. 


FMCD: 


FMC: 


ZM F 


MASK 


JU FETCH 


RMI 
JU $ 


FCOS 


AX, 0 


SFI 


FSIN 


FABS 


FEXT 


MR FMC:+ 3, AX 


JC AX, GEZ, FGET 


JU FACMP 


JU FGET 


LOC 


LOC 


FTBLE + FMCS 


FMC 


FTBLE + FMCSD 


FMGD 
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CHAPTER 6 


OPERATING INSTRUCTIONS AND SYSTEM GENERATION 


6.1 Using the Package as Supplied: 


A Command Equate Tape is supplied for the package as delivered. This 
tape is a source tape and contains the definitions of the floating point 
commands as well as the entries for S$SFI, the conversion routines, and the 

names and entry points for the various flags and routines in $SFI that can 
be used in non-interpretive mode or by user coded extended functions. A 
commented listing of this tape appears in Appendix B of this document. The 


tape supplied with the package is not commented. 


This tape is intended for use during assembly of the user programs 
which reference all or part of the Floating Point System. The user need 
not define the names or locations appearing on this tape when he generates 
the source tape for the program he has written. Assuming the source tape 
is«generated, the following procedure will assemble the program for use 
with SSFI: 

1) Load the Assembler and select pass 1. 

2) Position the Command Equate Tape in the reader. 


3) Press start. The assembler will come to a halt after reading the 


Command Equate Tape. 
4) Position the source tape for the program in the reader. 
5) Press continue. The assembler will now complete pass 1. 


6) Run pass 2 and pass 3 as usual, using only the source tape for the 
program. The Command Equate Tape need not be used during these 


two passes. 
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Steps 1-6 apply to each separate assembly of a user program. When 
the program(s) are assembled and are to be loaded, the user must be sure 
that SSFI is loaded before anything else, because any extended functions 
being used (including user-generated functions) overlay portions of FTBLE 


in SSFI. 


With the object tapes for $SFI, the conversion routines, and the 
extended routines, as delivered, it is possible to create one of the four 
memory maps depicted below. The user may generate his own system suiting 


his needs by following the procedures outlined in the next section. 
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POSSIBLE MEMORY MAPS 
WITH OBJECTS AS DELIVERED 
SSFI Only SSFI & Conversion SSFI & Extended 
Only Only Entire System 


7777. ‘| Utility 7777 Utility 7777. ‘| ~Utility 7777 Utility 
‘ Loaders Loaders Loaders Loaders 


i f 

4 

|  $SFI ; 

6120 | 

6117 | 

User i 

Area I 

771 | 
4770 


User 


Area 


User 


Area II 


0000 0000 0000 


0000 
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6.2 User Generated Systems; 


6.2.1 General 


The source tapes supplied for the basic package, conversion rou- 
tines, and extended package are organized so that the components of 
any Floating Point System the user may select for his problem will 
assemble such that the object tape generated will load as high as 
possible in memory. Since any configuration of routines which could 
comprise a Floating Point System must contain $SFI as a component, 
the procedure outlined for generating the system assumes that $SFI 
is located at the highest possible place in memory, and all other 


(if any) components will be located at successively lower locations. 


6.3 Assembling SSFI: 


Should the user desire to have S$SFI load at some other place in 
core, the package must be re-assembled. Since S$SFI will be assembled 
so that it occupies the highest location possible, the user need only 
supply that address. This is done by creating a separate source tape 


which has the following two assembly instructions: 


SSEND = XXXXX 


EOT 
where XXXXX is the last address that the user wishes SSFI to occupy. 


To assemble $SFI, one follows the usual assembly procedure, ex- 
cept that the short source tape constructed above must be read in before 


the $SFI source tape at each pass. 
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Note: Do not use the Command Equate Tape for this assembly. 


After completion of this assembly, the "SSFI=" statement on 


the Command Equate Tape should be updated (using the Source Text 


Editor) so that it States the new location at which SSFI begins. 


This location may be read off the symbol table printed out at the 


beginning of pass 3 of the assembly of SSFI. 


Generating a Floating Point System: 


If a configuration not obtainable with the object tapes supplied 


is desired, the user may, by editing the source tapes supplied (using 


the Source Text Editor), create a new source tape consisting of the 


routines he needs for his purposes. 


Since S$SFI is the integral component of any floating point con- 


figuration, it is assumed here that either the object tape supplied 


for SSFI 


sociated 


The 
arranged 
with the 


an index 


routines 


extended 


or a relocated version assembled as in 6.3 (and the as- 


revised Command Equate Tape) will be used. 


source tapes for the conversion and extended routines are 
so that there is one routine per block on the tape starting 
second block. (The first block contains comments comprising 


to the rest of the tape.) This makes it easy to extract the 


desired using the Source Text Editor. The routines on the 


package source tape are completely independent. However, the 


conversion source tape consists of three routines - the first (physically 


on the tape) of which is shared by both @SFC and @SCF. This means if 


the user 


desires a conversion routine in only one direction, he must 
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also include the shared routine (@FXC). The routines @SCF and @SFC 


are independent of each other. 


The procedure for extracting the desired routines from these two 


tapes is outlined below: 


1) Load and start the Source Text Editor 


2) Punch a single block at the beginning of the source tape 


being created which says: 
SSEND = S$SFI 


($SFI-1 will be the last location into which the resulting 
object tape will load.) 


3) If either or both of the conversion routines are desired, 
position the conversion routine source tape in the reader 
and skip the first (index) block. Copy the second block 
(@FXC), which is the shared routine; then copy the block(s) 
containing the desired conversion routine(s). If neither 


conversion is required, this step may be skipped entirely. 
4) If one or more of the Extended Package routines is required, 


position the source tape for the extended routines in the 
reader and skip the first (index) block. Then copy the 
blocks corresponding to the routines required. If no ex- 
tended package routines are required, this step may be 


skipped entirely. 
5) Punch a final block consisting of the single statement: 


END 
6) Load the Basic Assembler, select pass l. 


7) Position the Command Equate Tape which corresponds to the 


object of $SFI to be used in the reader. 


8) Press start. The assembler will come to a halt after reading 


the Command Equate Tape. 


9) 


10) 


11) 


12) 


13). 


Note: 
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Position the source tape generated in steps 1-5 in the 


reader. 
Press continue. The assembler will now complete pass l. 


Run pass 2 and pass 3 as usual using only the source tape 
created in steps 1-5. The Command Equate Tape need not be 
used for these two passes. (Even if no assembly listing is 
required, pass 3 should at least be started so that the 
Symbol Table is obtained. 


Check the Symbol Table entries for @FCST, @SFC and @SCF 
against the values appearing on the Command Equate Tape 
used for the assembly and, if they are not the same, create 
a new Command Equate Tape which has the values for @FCST, 


@SFC, and @SCF as they appear in the Symbol Table. 


Assemble any user programs using this new Command Equate 


Tape as explained in 6.1. 


1) The lowest location occupied by the Floating Point 
' System as created above will be the value appearing 
for SSEND on the Symbol Table printed as a result of 
step 11) above. 


2) The object tape created by this procedure must be 
loaded after $SFI as explained in 6.1. 


Code (octal) 


01 
02 
03 
04 
05 
06 
07 
10 


Definitions: 


[D] ~optional selection of deferred addressing 


00 
[101] 
[102] 
[103] 
[104] 
[105] 
[106] 
[107] 
[110] 

14 

15 

16 

17 

20 

21 

22 


Y ~ address of floating operand 
~ address of location containing address - 1 of floating operand 


~ address of another floating command 


APPENDIX A 


Command Summary = Basic 


~ address of index value 
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I ~ index value of source or destination at address Y 


A ~pseudo-accumulator (FAC) 


X - pseudo-index register 


F ~ floating value of source or destination at effective address formed 


from Y. 


FEXT 


FLDA [D] 
FSTA [D] 
FADD [D] 
FSUB [D] 
FMPY [D] 
FDIV [D] 
FADM [D] 
FSBM [D] 
 ‘FABS 
FASQ 
FNOR 
FNEG 
FJMP 
FJAP 
FJAZ 


Basic Commands 


ee we Ke KKK OK 


Operation 


exit 
FA 
A->F 


ATF ~oA 


A-F —bA 
AXF—> A 
A/F —-PA 
A+|F|>A 
A-|F]—-»A 


[ala 


A2—>A | 


normalized A~PA 


-A—PA 


jump to Y 


jump to Y if A>O 
jump toYifA#0 


Flags 


none 
none 

FXFLG 

FXFLG 

FXFLG 

FXFLG 

FXFLG, FDFLG 
FXFLG 

FXFLG 

none 

FXFLG 

FXFLG 

none 

none 

none 


none 


Registers 
none 
FAC , FTM 
FAC 
FAC, FIM 


- FAC, FIM 


FAC, FIM 
FAC, FTM 
FAC, FTM 
FAC, FTM 
FAC 
FAC, FIM 
FAC 
FAC 
none 
none 


none 


Code (octal) 


23 
24 
25 
26 
27 
30 


Code (octal) 


(*) 


(+) 


(1) 


31 
32 
33 
34 
35 
36 


If 
be 


If 
be 


If 
be 
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Basic Commands Operation Flags Registers 
FJAN Y jump to Y if AZO none none 
FJEV Y jump to Y if FXFLG set#0 FXFLG (set to zero) none 
FJDC Y jump to Y if FDFLG set£0 FDFLG (set to zero) none 
FJIX Y X+1-®X, jump to Y if ¥/0_ none FINDX 
FLDX Y I—>Xx none FINDX 
FSTX Y X ~el none none 
Command Summary - Extended Functions 
Extended Command Operation Flags Registers 
FSIN SIN (FAC)—® FAC FXFLG FAC, FTM 
FCOS COS (FAC) —®* FAC FXFLG FAC, FTM 
FATN TAN-1 (FAC) —®FAC FXFLG,FDFLG FAC, FTM 
* 
FLNE LOG, (|FAC|)—» FAC rxFic**? FAC, FTM 
FEXP eFAC_» Fac FXFLG, FoFLG‘)) pac, Frm 
FSQT | FAC] =» FAC none (+) FAC 


input argument is 


set non-zero (see 


input argument is 


set non-zero (see 


input argument is 


set non-zero (see 


write-up). 


write-up). 


write-up). 


negative, FNLNF internal to the FPLNE routine will 


negative, FSFLG internal to the FPSQT routine will 


too large, FEXOF internal to the FPEXP routine will 
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01 3$FCQ - COMMENTED pee ee MeN 
eve 3 74-43-4021 far ree 
B03 3 COMMAND EQUATE TAPE 

O04 $SFI=6120 3 ENTRYs INTERPRETEX 

05 @FCST=5765 3 CHARACTER SET TABLE 

wv6 @ SFC=5377 3 ENTRY, FLING. TO CHAK. 
oT @SCF=4771 3 ENTRY, CHAK. TO FLING. 
W008 ¥GET=$SFI+2 3 EN1kYs EXTERNAL RETURN 
boy FM ASK=FGET+12 3 MASK FOK DEFERRED ODE 
010 FPUNT=FMASK+13 3 ILLEGAL COMMAND ROUTINE 
11 FTBLE=FPUNT+14 3 TABLE OF ENTRIES 

le FTCMP=FTBLE+ 111 3 ENTRYs Fim NEGATION 

013 FTCMA=FTCMP+4 3 ENTRY, -CAXsAY) TO FTd 
14 FTSAX=FTCMA+5S 3 ENTRYs CAKsAY) TO FIM 
15 FIMHI=FTSAX+1 3 HIGH ORDER FT 

016 FTMLO=FTMHI +2 3 LOW ORDER FT 

B17 FACMP=FTMLO+2 3 ENTRYs FAC NEGATION 

018 FACMA=FACMP+4 3 ENTRYs -(AKsAY) TO FAC 
wig FASAK= FACMAt5 3 ENTRY, C(AKsAY) TO FAC 
vey F ACHI=FASAK+1 3 HIGH OKDEH FAC 

g21 FACLO=FACHI +2 3 LOW ORDER FAC 

22 FNORM=FACLO+2 3 ENTRYs NORMALIZATION 
923 FXFLG=FNOHM+3 3 EXPONENT OVERFLOW FLAG 
B24 F ACXP= FXFL G+ 50 3 FAC EXPONENT 

B25 FINDA= FACKP+20 3 PSEUDO- INDEX 

v26 F DAD= FINDA+5 3 ENTRY» DBLE. PREC. ADD 
027 FDMPY=FDAD+ 31. 3 ENTRYs DBLE. PREC. MULT. 
Bau FDDI VsFDMPY+143 3 ENTRYs DBLE. PHECe DIV. 
Bey FDFLG=FDDIV+3 3 DIVIDE CHECK FLAG 

630 F SDVD=FDFLG+ 121 3 ENTRYs SNGL. PREC. DIV. 
031 FLO DV=FSDVD+14 $ Le O« DIVIDEND FOK ABOVE 
032 FETCH= FLO DV+27 3 ENTRY. FLING ARG TO FT 
633 F SPL T=FETCH+ 16 3 SPLIT ARG TO FAC. 

034 F IMXP= F SPLT+6 3 FT EXPONENT 

035 FAKGD= FTMAP+24 3 ENTRYs IMMED FETCH TO Ax 
036 FPUSH= FARGD+1 3 PUSH DOWN POINTER 

037 FARSN= FPUSH+2 3 ENTHYs ARITH RGHT SHF1 
238 F ARSC=FARSN+6 3 SHIFT COUNTER FOR ABOVE 
639 FCMAX=FAKSC+7 3 ENTRY, MAK. NU. TO FAC 
vay FOFAC=FCMAK+21 3 ENTRY, ZEKO TO FAC 

Q41 FPSTA=FQFAC+ 32 © 3 STOKE FAC 

re FINAC=FPSTA+34 3 PACK FAC 

B43 FLI ST=FINAC+373 3 PUSH DOWN LIST 

B44 FEXT=@ 3 EXIT COMMAND 

045 FLDA=1 3 LOAD FAC CoMmmAND 

046 FLDAD=101 3 LOAD FAC DEFERRED 

47 FsTAa=2 3 STORE FAC 

048 F STAD= 102 $ STORE FAC DEFERRED 

gay FADD= 3 3 FLOAIING ADD 

w5e FADDD=103 3 FLOATING ADD DEFERRED 
051 FSUB=4 3 FLOATING SUBTRACT 

use F SUBD=104 3 FLOATING SUB. DEFERRED 
53 FMPY= 5 3 FLOATING MULTIPLY 

O54 FMPY D= 105 3 FLOATING MUL1+ DEFEKHED 
055 FDI V=6 3 FLOATING DIVIDE 

056 FDI VD= 186 3 FLOATING DIVIDE DEFERRED 
057 FADM=7 3 FLOATING ADD MAGNI TUDE 
058 FADMD= 107 3 FLING ADD MAG DEFERRED 
ws9 F SBM= 10 3 FLOATING SUBe MAGNITUDE 
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660 FSBML=A11¢ 3 FLING SUB MAG DEFERREU 
061 FlAKN=11 3 IRACE UN 
W62 FinkF=le > IHACK OFF 
063 FSET=13 > SE! ERROR inAP 
v64 FABS= 14 > ABSIJLUTE MAGNI TUDE 
265 rAsw=15 > SQUARE 
066 FNOn= 16 > NIAMALIZCE 
067 FN&G= 17 3 NEGATE 
068 FJMr=2eV > UNCOUNDITIVNAL JUMP 
069 FJAP=21 > JUMP IF FAC > JR = ¥ 
870 FUAZ=2e2 > JUMP IF FAC = wb 
071 FUAN=23 > JUMP IF FAC < ¥ 
W772 FULVE2S3 > JUMP IF FAFLG WOl wW 
873 FJLC=25 3 JUMP IF FDFLG NU oY 
O74 FJIA=26 > BUMr FINDAs JUMr IF NOI 
075 FLDA=27 > LUAD PSEUDJIU-INDEA 
v76 FSTA= 30 5 O10nE POSLULJU-INDEA 
O77 FSIN=31 > oINE 
VIS FCOS= 32 > CUSINE 
O79 FAIN=33 > AnCIANGENT 
ron) FLNE=34 3 NAIURAL LOGAnI in% 
Vd 1 FEAP=35 > BAPOUNENITIAL 
Wdée FSW l= 36 > SQUARE KOJI 
O63 Bl > ASSEMBLER HALT 
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APPENDIX C 


Source Tape Organization 


NAME: SSFI 


3SSFI 


374-43-401L 

; IDENTIFICATION 
$SEND=S$SEND-1467 
LOC S$SEND 


block mark 


comprises several blocks 


‘block mark 


block mark 


3 @FXV 
_ 374-43-410L 
| p IDENTIFICATION 
' s INDEX (BLOCK NO. 
' VS. ROUTINE) 


block mark 


 $@FCST 

-374-43-411L 

; IDENTIFICATION 
$SEND=$SEND- 213 
LOC S$SEND 
@FCST=. 


block mark 


;@SFC 
374-43-412L 

; LDENTIFICATION 
SSEND=$ SEND- 306 
LOC $SEND 
@SFC=. 


block mark 


eo 


| 
'3SCF | 
1374-43-413L 
';IDENTIFICATION | 
_ $SEND=SSEND- 406 
LOC SSEND 
@SCF=. 


block mark 


eee oe et 


END 


' ‘i 
H er Se ei, ss Se te 


block mark 
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\ 
| 

‘ index (block 0) 

\ 

i 

common module (block 1) 

1 

/ 

} floating to character (block 2) 
y, 


2 


character to floating (block 3) 


NAME;  $FEF 


( 


Rese seria sor -aaneneegn oon say pean 


;SFEF 

374-4 3-420L 

; IDENTIFICATION 
;INDEX (BLOCK NO. 
' 3VS. ROUTINE) 


| ou 


ae eR OR A AC +t 


, 


block mark 


: 3FPSIN | 
374-43-421L 
_ 3 IDENTIFICATION | 
 §SEND=$$END-173 | 
LOC $$END 
FPSIN: | 


block mark 
poorer 
_ 3FPATN 
' 374-43-4221L 
SSEND=$ SEND-146 
LOC SSEND 
FPATN: 


block mark 
ee TEE 
' 3FPLNE 
: 374~43-423L 
| 3 IDENTIFICATION 
| S$SEND=$SEND~164 
/ LOC $$END 
| FPLNE: 
bo ee 
block mark 
Pinatas Bay 
_ 3 FPEXP 
' $74-43-424L 
' 3 IDENTIFICATION 
$SEND=$ $END-201 
LOC $SEND 
FPEXP: 


block mark 


i 
+ 
? 


3FPSQT | 
374-43-425L | 
; IDENTIFICATION 

SSEND=$SEND~127 
LOC $S$END ; 
FPSQT: © 


oe me 


block mark 


[om 


i 


Fee ee et 


block mark 


a 


oe ad 
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index (block 0) 


Sine/Cosine (block 1) 


arctangent (block 2) 


natural logarithm (block 3) 


exponential (block 4) 


square root (block 5) 
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APPENDIX D 


“FCG - Floating Point Constant Generator 


“FCG — Floating Point Constant Generator 


“FCG is a utility routine which is provided should the user wish to 
use floating point constants whose octal equivalences are unknown. With 
“FCG, the user can type in a floating point decimal number and receive 


the equivalent internal floating point representation. 


#FCG occupies locations 0-2660 inclusive. 


Operating Instructions 


1. Load “FCG by means of ZALD. 


4. Press START. 
5. FCG responds with a carriage return, line feed. 


6. Type a string of up to 1346 characters terminated with an equal 
sign (=). The character string should be in the format described 
in section 3.3, where the delimiter is an = rather than a carriage 
return. Typing a back arrow at any point causes the first previous 
non-back arrow to be ignored. Typing rubout at any point causes 
“FCG to type a carriage return, line feed, question mark (?) and 
returns to step 5. Typing more than 13 characters before typing 


an equal sign has the same effect as typing rubout. 


7. When the user terminates the character string with the equal 
sign, “FCG responds by typing the 2 word floating point equivalent 


(in octal) and returns to step 5. 


8. If the character string did not conform to the format specified in 
section 3.3, the message SCAN ERROR is typed and “FCG returns to 
step 5. 
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If the character string resulted in a number whose magnitude was 
outside the range 1.469369E-39 to 1.701411E+38 or if the character 
string contained more than 1059 mantissa digits, the message 


ANSWER OUT OF RANGE is typed and “FCG returns to step 5. 


E~-1 
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APPENDIX E 


FPSET —- Error Trap Routine 


Introduction: 


A series of floating point calculations on an unknown data base can 
generate errors, such as results which exceed the capacity of the machine 
or dividing by 9, etc. In order to facilitate the localization of the oc- 
currence of such errors, FPSET is provided and serves as an error trap 
routine. When an error specified by the user is detected, FPSET will 
interrupt the operation of the interpreter and give control to a user 
supplied error routine. FPSET supplies the user error routine with the 
following information, allowing the user to pinpoint the step in his cal- 
culations at which the error occurred: 

AX = recursion level at which the command at the address in AY was 

executed. 


AY = address of command executed immediately previous to detecting 


the error. 

TRP = error number indicating which flag in the user supplied error 
list was set non-zero. (TRP = position of address of error 
flag in user supplied table (see usage)). 

The recursive capability of the interpreter somewhat complicates cer- 
tain usages of FPSET and, for this reason, three modes of operation of 
FPSET are allowed: "On", "Off", and "Partially On". The latter mode 
allows FPSET to keep track of commands and recursion levels without ex- 
amining any error flags. The utility of this mode is described in the ex- 


amples at the end of this appendix. 
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Usage: 


FPSET is controlled by the use of the FSET command in the sequence of 
floating commands being executed by the interpreter. There are three modes 


of operation of FPSET: 1) ON, 2) OFF, and 3) PARTIALLY ON. 
1) To turn FPSET ON, the command is 


FSET A 


where A is the address of a table with the following format: 


A: WRD ERR ;USER ERROR ROUTINE ENTRY 
WRD FLGL 3ADDRESSES OF SYSTEM 
WRD FLG2 ;FLAGS TO BE CHECKED... 
WRD -1 3;END OF TABLE SIGNAL 


When this FSET command is encountered with a positive non-zero value for 
A, FPSET will examine the state of every flag listed in the table at address A 
after every command executed by the interpreter from the point of the FSET A 
command onward. Whenever a flag whose address is in the user list has become 
non-zero (indicating an error), FPSET zeros the flag, and then gives control 
to the user error routine at the address specified in the first word of the 
table at A. The information supplied to the user error routine is as stated 


in th introduction. 


The user error routine may use $SFI, but any additional errors which might 
occur will not be checked by FPSET, and any FSET commands in the command se- 
quence will be ignored. If the user wishes to call $SFI in his error routine, 
it is up to him to save and restore the states of the interpreter system flags 


and the floating accumulator (FAC) before and after such SSFI use. 


FLOATING POINT MANUAL 
71-44-001 
2) To turn FPSET "off", the command is 


FSET 0 
This completely disconnects FPSET from the interpreter. 
3) To turn FPSET "partially on", the command is 


FSET N 


where N is any negative number. 


In this mode, FPSET will keep track of the current command address and 
recursion level but will not examine any flags. If FPSET is at some later 
time turned "on" and discovers a flag set non-zero, the level and the com- 


mand address will be correct within certain limitations (see Notes). 


This mode is useful when the user does not wish to enter his error 
routine for errors which occur during execution of a section of his command 
sequence. For example, the command sequence may contain an FJEV or similar 
test for conditions known to the user, and with FPSET "on", these conditions 
could be altered (cleared) if the corresponding flags are in the user error 
list at A. In this case, an FSET N (where NCO), issued before entering this 
section, and an FSET A (ADO, A=address of table), issued after completion of 


this section will allow FPSET to retain the necessary information should 


other errors occur and allow the section itself to operate properly. 


User Error Routine: 


Basically, the user error routine may do anything. However, the user must 
remember that his error routine is considered as an extension of the inter- 
preter. At the completion of the error routine, control should be given back 


to the interpreter via a JU FGET or similar return. 


Register AX is used as an argument upon return (via JU FGET) and can turn 
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FPSET "on", "off", or "partially on" according to AXS0O, AX=0, or AX<0 


respectively. If AX>0, it must be the address of an error table as de- 


scribed above (it need not necessarily be the same one as before). 


Notes: 


1) 


2) 


3) 


It is generally the case that an error flag is set by the command 
immediately preceding the detection of the flag non-zero. In the 
case where FPSET was not "on" at that moment, but was turned on 
later and found the flag non-zero, FPSET will report that it does 
not know which command caused the error by giving an AY value 
which points to the FSET "on" command or by AY = -l. The dif- 


ference in meaning of the two AY values is as follows: 


a) AY = address of FSET "on" command if the flag was non-zero 


at the time the FSET "on" command was encountered. 


b) AY = -1 if all flags were zero when the FSET "on" was en- 
countered, but a flag was set non-zero later at a point 
which indicated that the command which caused the error 
was at a recursion level one less than the level at which 
the error was detected. This situation is avoided if 


FPSET is partially on throughout until it is turned on. 


If the user wishes to restart his entire program or in any other way 
wishes to use the interpreter without reloading it, he should make 
sure that FGET+l is initialized to FARGD and FPUSH is initialized 

to FLIST-1. 


A FEXT command does not affect the mode of operation of FPSET, i.e. 
upon re-entering $SFI, FPSET will operate as per the last FSET com- 


mand encountered before the FEXT. 


Operating Instructions: 


oi 


If a user generated floating point system is being used (as per 
section 6.4), the FPSET program needs to be reassembled. To do 


this, make a short source tape with the following two instructions: 


SSEND = XXXXX 


EOT 
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where XXXXX is the lowest location occupied by the user generated 
system. Then assemble FPSET as explained in section 6.1 of this 
document, making sure that the short source tape above is read in 
before the FPSET source at the beginning of each pass (steps 4 & 6 
of 6.1) and the Command Equate Tape is read in at the beginning 
of Pass l. 


Load FPSET after SSFI. 


Start user program which has FSET commands in the usual way. 


Examp 


les 


1. Typical usage of FSET "on" and "off" 


U 


ser Main Code 


JU $SFI 
FSET A 
FD 

FSET 0 
FEXT 

WRD ERR 
WRD FXFLG 
WRD FDFLG 
WRD -1 


3Turn FPSET on User Error Routine 
causes divide check, FPSET gives contiro O ERR: . (TRP = 000002 
This returns control to interpreteh . indicating second 
flag in table at A 


to execute next command : was set) 


Errors in here . 
are not detected MRI A,AX 


by FPSET ;Tlurn FPSET back on 
y 


JU FGET 


;Address of user error routine 
;Address of exp over/underflow flag 
;Address of divide check flag 

3;End of table 


| 


T00-77-TZ 
TVONVN INIOd ONILVOTS 


Example 


2. Use of FSET "Partially On" 


User Main.Code — User Coded Extended Function 


JU $SFI | ae FPUFN: JU  $SFI 


;Turn FPSET on 
(Detects FDFLG error 
immediately and traps 
to user error routine) 
;Reason for FSET partially on 
3sError (divide check) occurs here sUser supplied table 


;User defined function 


In this example, because FPSET was only partially on when the error actually occurred, the 
error trap will indicate that the erring command was FSET. It will, however, indicate that the 


divide check flag was on (TRP = 000001 since FDFLG is first in table A) and that the error oc- 


curred in recursion level 2 since FPSET was turned on in the user function. 


T00-79-TZ 
TVWONVN INIOd DONILVOTI 


Examples 


3. Another usage of FPSET "Partially On" 


User Main Code 


User Coded Extended Function 


JU SSFI 


JU SSFI 


3Turn FPSET partially on 


;lurn FPSET on 
(no errors yet) 


3;User defined function 


In this example, FPSET indicates an error on 


Input argument cause 
user coded function 
to set error flag 


recursion level 1, and that the command 


ZM P1,FUFLG 


causing tie error was FUFN (i.e. AY will have 
address of FUFN command). This is as it should 


be since the arguments given to the user function at 


error is detected here 


the FUFN caused the function to set an error con- by FPSET 


dition. 
Note: If the FSET -1 had not been issued in the Main 
Code, FPSET would have indicated the AY= -1 condition. 


It would, however, indicate the correct flag and the 


correct level (i.e. level 1). WRD FDFLG 
WRD FUFLG 

: 

WRD -1 


T00-77-TZ 
TVONVW INTOd DNILVOTA 


8-4 
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APPENDIX F 


Trace Routine 


The floating point trace is a debugging aid which prints 
the value of pertinent variables in $SFL and the user's program 
before the execution of each floating point pseudo command. The 


variables printed are: 


A. current level of $SFI 

B. address of the instruction to be executed 
C. code for the instruction to be executed 
D. FINDX (floating point index) 

E. FDFLG (divide check flag) 

F. FXFLG (exponent overflow/underflow flag) 
G. FAC (floating point pseudo accumulator) 
H. effective address of argument, if any 


I. value of argument, if any 


The user specifies which of the variables are to be printed 
and the maximum level for which he wants the information printed. 
This is done through the floating point pseudo commands FTRN 


and FTRF. 


To turn the trace on and specify which of the nine variables 
are to be printed, the pseudo command is: 
FIRN X 
where bits 0-8 of the integer X correspond to the variables A through 


I above. For each bit that is on (=1) the corresponding variable 
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will be printed before the execution of each floating point command. 
The FTRN command sets maximum recursion level to be traced tc 7, 

turns the trace on and prints a heading (A-I), telling which variables 
are to be printed. The "trace on" causes the specified variables to 


be printed on one line before each instruction is executed. 


The printed value of variables H (argument effective address) 
and I (argument) need further explanation. If tte command to be 
executed has no argument, columns H and I will be blank. If the 
argument is floating point, I is printed as a floating point deci- 
mal number, otherwise it is octal. If the command is FTRN, FTIRF or 
a JUMP command, then H is the address+l of the command and I is the 
contents of H. In the case of the commands FLDX Y, or FSTX Y, H is 
the address Y and I is the contents of Y. For user coded extended 
functions, columns H and I will be blank unless the function has a 


floating point argument. 


To turn the trace off beyond a certain level, the pseudo command 
is: 
FTRF X 
where the integer value X specifies the maximum recursion level (1-7) 
for which the specified variables are to be printed. If X is less 
than or equal to 0, the trace is disabled and no variables will be 


printed from then on until another FTRN X command is executed. 


FLOATING POINT MANUAL 
71-44-001 


Notes: 


1) When the trace has been turned on, certain locations in 
SSFI are changed. SSFI is restored to its original state 
only after the trace is completely disabled by an FTRF 0 
command. Therefore, to restart the user program or use 
SSFI without reloading when the trace has been on, the user 


should make sure that: 


a) FGET+l1 is initialized to FARGD 
b) FSPLT is initialized to 11 0000 06 
c) FSPLT+1l is initialized to FTMHI 


d) FPSTAt1l is initialized to FARGD 


2) The trace program cannot run at the same time as FPSET, 
and so should be assembled over the FPSET program. Therefore, 
if the trace is on and an FSET command is encountered, the 
FSET command is considered illegal. That is, the interpreter 
haits with the address of the iliegai command in the MB reg- 


ister on the front panel. (See section 2.2.5). 


Operating Instructions: 


1) If a user generated floating point system is being used 

(see section 6.4), the trace program (FPTRC) should be re- 

assembled. Since FPSET and FPTRC are mutually exclusive, 

FPTRC is assembled over the FPSET program. To do this, make 

a short source tape with the following two instructions: 
SSEND = XXXXX 


EOT 


F-4 
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where XXXXX is the lowest location occupied by the user 
generated system (excluding FPSET). Then assemble FPTRC 
as explained in 6.1. Make sure that the short source tape 
above is read in before the FPTRC: source at the beginning 
of each pass (steps 4 & 6 of 6.1). Read the Command 


_Equate Tape in at the beginning of Pass l. 


2) Load FPTRC after SSFI, 


3) Load and start user program with FTRN and FTRF commands 


in the usual way. 


061 
W082 
083 
O84 
WS 
006 
O07 
B08 


089 


012 


ae os SE 


OBO08 
00001 
68002 
00083 
G0G04 
06085 
00GK6 
BO087 
06016 
08011 
O801e2 . 
66013 
66014 
Q90B15 
00616 
QOB17 
YBV20 
06021 
00622 
0O023 
BOG24 
00025 
06826 
880627 
66030 
68831 
606832 
60033 
66034 


28835 


a7 as 


06 8180 
606126 
BB B8800 
OBO757 
0B B60 
090832 
BO OWOH 
0868633 
00 0080 
00080635 
04 BBB 
068835 
00 BBB 
YOBORE 
66 O8B08 
080367 
BO BHO 
908033 
64 6806 
060035 
66 8008 
04689037 
64 8000 
0400008 
08 BOBO 
02 0108 
177775 
B50000 
800263 
B4ABB88 
OB020e2 
OBBBBB 
BBBVBO 


5 EXAMPLE 1 
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sTHIS EXAMPLE USES TWO FTN COMMANDS-~ 
$THE 15ST FTRN SELECTS VARIABLES A»sBs Cs Ds Fs GsHs I 


370 BE PRINTED. 


THIS IS IN EFFECT UNTIL THE 
32ND FTRN IS EXECUTED. 
SARE THEN CHANGED TO As Bs Cs Eo Fs GoHe 


THE SELECTED VAHI ABLES 
THE 


- SFTRF @ COMPLETELY DISABLES THE TRACEe 


03 


11 


27 


Ql 


O35 


Be 


26 


11 


01 


66 


i2 


le 


OB 
06 


$SFI 


7357 


FLDA W 


HLT 
-3 
98OB0s 203 


LWP ATATAAA. OMO 
Ys: imp 49008»292 


es) 


SPRINT As Bs Cs Ds Fs Go Hs I 


3 Y=A*Y 


> DONE LOOP 3 TIMES? 


SPRINT As Bs Cs Es Fs GH 


3s TURN TRACE OFF 


3LOOP COUNT 
3508 


° fA 
7GOe VY 


3 O08 


ee ee ee en ee ee ee od Me 


r 


a a ae 


40801 
OBB806 
WdBW2 
WHble 


D 


YVOUOHO 
177775 
177775 
177775 
177775 
177776 
177776 
177776 
177776 
177777 
177777 
177777 
177777 
OYWbUO 


E F 


®e ese 
Gece 


‘ry 


SSes Seesqegegcgecseck 


F-6 


G 


+20 35622TE- 39 
+20 356227E- 3y 
+5eOUBUUOE+ OU 
+1-G00080E+O1 
+1-¢00Q000E+G1 
+1-000000E+01 
+52 OBGOW0UE+ 00 
+52e¢B00000Et+01 
+5-¢O60000UEtTO1 
+5-¢OUVBUOOE+O1 
+5eWUVOWUE+OU 
+206 SOUOUBE+ 2 
+220 500000E+Ge 
+20 SUYUHWUEtT Oe 


G 


tee SUUWUUE+U2 
+5-OGBHBWUEt OO 
+2eW00HW0UUE-bKe 
+2-00Q000E- 2 


H 


O0V3e 
00033 
0Y0W35 
00835 
080815 
06033 
GOB35 
BUB35 
BOW15 
60033 
00035 


BBYB35. 


BBO15 
00017 
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I 


177775 
+5eGOWUUYE+OU 
+2eBUKOUUE+O 
+2eW0UO0UUEt+ OU 
8060606 
+5eBBGOOUWE+HO 
+1-GO0Q0000E+GO1 
+1-Q@0000UEt01 
OUW006 
+5eWUVUBWUVUETO 
#5-OW0BO0UEtO1 
+5e¢UW0W00UE+01 
YVOBLYE 
000367 


YVWWOB 
68001 
Bbdbe 
BO063 
OBBQB4 
BBWOBS 
OBUO6 
OBO07 
BYH16 
60011 
w00le 
06013 
00014 
00615 
00616 
BYB17 
BWVD2O 
00021 
OO0e2e2 
00023 
40624 
00625 
00826 
OB027 
00030 
06831 
WOO32 
98033 
00034 
00035 
60036 
00037 
BYO048 
00041 
06642 
00043 
06044 
ZBB45 
QOBB4A6 
668647 
08050 
00051 
06160 
00161 
dvide 
00163 


B62e¢d2 


Wd 01GB 
046120 
00 OBB 
O40777 
08 BY0H 
0OBB0G2 
00 BOW 
VUBB43 
QB 0601 
OBBB44 
06 6800 
080646 
YH YOO 
060046 
00 WYO 
000010 
00 8000 
QOBBBO 
BB 89608 
be 0100 
06 8109 
007057 
il 6860 
YUB035 
0B W1BY 
0¥61le0 
04 8008 
BHB058 
04 Jd 
VUOBBOG 
0B 6001 
Q6O0UW45 
00 8006 
00 0108 
O¥61ee2 
177776 
BOGBT7 
000077 
000800 
BUUBOO 
076406 
000211 


862000 
008266 
BS50080 
BBB204 


VOdB24 


00 
03 


3 BAAMPLE 2 
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sTHIS EXAMPLE PRINTS ALL 9 VARIABLES 


3FUn RECURSION 


LEVELS 1 AND ae 


3 THE USEK EATENDED FUNCTION Is AT 


sHECURSION LEVEL Ge 
JU $5FI 


FITAN 777 


FTRF 2 


FLDX Xx 


FLULAD X+1 


FRPUFNS 


FLDA Z 


AHG: FDIV @ 


FSTAD Ate 
FEAT 
JU 

WHD =< 
Wnb 
WRD 
WixD 


WnhD 764005211 


100 
62000» 206 


LOC 
WH 
WKD 360W0s204 
FUFN= 37 
LUC 
WrD 
EN D 


PYBLE+FUFN 
FPUFN 


SPrRINY As Bs Cs Da Ko Fs Goris I 


3MAX TRACE LEVEL=2 


SFPEICH ARG DEFERKED 


SUSENK LATENDED FUNCTION 


> DINE? 


SYESs TRACE OFF 


+ 50GBe¢O/Y 


sLO0OP COUNT 
3 FETCH ADHs 
3 STJRE ADRs 
3FUFN ARGe 


DEFERRED 
DEFERRED 


3300-0 
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BOBO 
66001 
68002 
YYW 3 
YWOWOA 
BYBBHS 
BVOBU6 
Bd007 
AWU1 
80011 
68012 
00013 
GYVG14 
80015 
68016 
00017 
6UB20 
80021 
82022 
Hbd23 
O06024 
00625 
A4026 
80827 
06030 
04031 
OBb3e2 
BBB33 
BBB34 
66035 
BdB36 
86837 


BBVA 
66041 
0OV42 
00043 
BBO44 
OBW45 
OBH46 
OBVB47 
Y4050 
606051 


40168 
90181 


60102 
00163 


06222 


66 0100 
806128 
0B B00 
000777 
00 OBOE 
BOGBO1 
06 BBBE 
0208043 
0) @B61 
JISBG4S4S 
02 BBBG 
0BBG46 
06 BHU 
OBOB46 
00 BOBO 
0BB018 
BO BBO 
BBBGBG 
48 O40 
Y2 0106 
08 0100 
667657 
11 2800 
080635 
00 8160 
606120 
06 BOO 
600058 
88 8080 
B6BBK06 
60 6061 


WOW 4 
68 BBB 
64 0108 
646122 
177776 
OOOSTT 
000077 
OOUBVWO 
BOBBLO 
076486 
060211 


662606 
VBO0206 
850008 
680204 


600824 


03 
1 
12 
27 
01 
G2 
37 
26 
12 
0G 
Qo 
03 


06 


03 


01 
86 
v2 


0@ 
63 


3 EXAMPLE 3 


sTHIS Is THE 


3 REACEPT 
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SAME AS EXAMPLE 2 
THE MAXIMUM RECURSION 


SLEVEL PRINTED IS LEVEL 1 DUE 


3TO THE FYTRF 1 COMMAND. 


JU $SSFI 
FTRN 777 
FTRF 1] 
FLoOn X 
FLDAD XA+1 
FSTA Y 
FUFN Y 
FJIA e-6 
FIRE @ 
PEAT 
FOM HLT 
FPUFN: JU FARGD 
RY AXs AXG+ 1 
JU $SFI 
FLDA Z 
AnGs FDIV @ 
FSTAD Ate 
FEAT 
JU FGET 
As WRD -2 
WHD 77 
WRD 77 
Ys WRD 656 
Zs WRD 76400;211 
LUC 1800 
WRD 620060> 206 
WRD 560003204 
FUFN= 37 
LOC FTBLE+FUFN 
WRD FPUFN 


EN D 


SPRINT As Bs Cs Ds Es Fs GsHol 


sMAX TRACE LEVEL=1 


3 FETCH ARG DEFERRED 


3USER EXTENDED FUNCTION 
3 DUNE? 


3YESs TRACE OFF 


3 3O0B-0/Y 


3LO0P COUNT 
3 FETCH ADRs 
$ STORE ADRs 
3 FUFN ARG. 


DEFERRED 
DEFERRED 


3568-0 


I> 


.8 


a) 


ROE by 
1 CATAB 
ROOVE 
CRYO 
COETA 
ARP TA 
BACT CA 
CEO D 
(Guancan re 
(SES 7, 


APA eoe 
(1 AO OG 


C 


eerLeS 
OAreFT 
AR TAY 
AReRO 
PRART 
CEPA 
ROLE Y 
PRAOD 
Pe eORT 
PCROE 


AOL? 
Pe] 


4 


CHAAR 
ECAC RAE 
177776 
ie Bt a aE 
177774 
PET EG 
177774 
Ni a a a 7 
TTT T 
LET LY 
CLP ONIN CRIA 


++ t+ bt te + + + 


F-10 


-_~ 
$4 


avs 
e 


REEL TRH BO 
S5620F The 49 
B562 Oo The BO 
oe CACEARKE G] 
Se MP MAAE RCI 
LeFRPCEARSE I 
Le PACA OAR+ 2 | 
LehECARER+G 1 
Le AGCOCERLE 
5. FAPKHACKE EY 
BSCE OMAA EET 


“ 
aA NG 
e 


Nfs 
e 


mn 


FLOATING POINT MANUAL 


H 


BAGG 
Ere aAS 
iam oad 
AO OLE 


PEO 
GREE 
COLE 


CECE 


71-44-001 


COORE | 
77776 
tS. FOAEMOEEOY 
FO CPOE OETE NE 


fPEALE 
tLe CARE ES 
tf EELR EEA | 


Wa] 


COBRA 


Sa 3 " ~~ aay 
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